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Effects of two types of shock topology, namely, small-scale shocklet and large-scale shock 
wave, on the statistics of temperature in compressible turbulence were investigated by 
numerical simulations. The shocklet and shock wave are caused by the purely solenoidal 
and primarily compressive modes of large-scale random forces, respectively. Hereafter, 
the corresponding two flows are abbreviated as SFT (solenoidal forced turbulence) and 
CFT (compressive forced turbulence), respectively. It shows that in SFT the tempera¬ 
ture spectrum follows the fc -5 / 3 power law, and the temperature field has the ” ramp-cliff” 
structures. By contrast, in CFT the temperature spectrum defers to the /c -2 power law, 
and the temperature field is dominated by the large-scale rarefaction and compression. 
The power-law exponents for the probability distribution function (p.d.f.) of large neg¬ 
ative dilatation are —2.5 in SFT and —3.5 in CFT, very close to those computed from 
a theoretical model. In CFT, the collapse of the p.d.f. for temperature increment to the 
same distribution indicates the saturation of scaling exponent at high order numbers. 
For the isentropic assumption of thermodynamic variables, it shows that the derivation 
in SFT grows with the turbulent Mach number (M t ), and for same M t , the variables 
in CFT are more anisentropic. The angle statistics of CFT shows that the temperature 
gradient is preferentially perpendicular to the anisotropic strain rate tensor. In detail, it 
tends to be parallel with the first eigenvector and be orthogonal with the second and third 
eigenvectors. By employing a ” coarse-graining” approach, we investigated the cascade of 
temperature. It shows that the temperature variance at large scales is increased by the 
viscous dissipation at small scales and the pressure-dilatation at moderate scales, but is 
decreased by the subgrid-scale (SGS) temperature flux, which preferentially transfers in 
the orientation where the temperature gradient anti-aligns with the SGS temperature- 
velocity coupling. The distributions of pressure-dilatation and its two components prove 
the fact that the negligible contribution of pressure-dilatation at small scales is due to the 
cancelations of high values between rarefaction and compression regions. The strongly 
positive skewness of the p.d.f.s of pressure-dilatation implies that the conversion from 
kinetic to internal energy through compression is more intense than the opposite process 
through rarefaction. Furthermore, it shows that in SFT the fluctuations of pressure- 
dilatation approximately follow the Zeman model (Zeman|[l991). 
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1. Introduction 


Since the earlier seminal work of Corrsin (Corrsin 1951), much effort has been devoted 


to studying the statistics of temperature fluctuations in turbulence, including astro¬ 


physics (Cattaneo et al. 2003), geophysics (Stevens 2005) and engineering (Hill 1976). 


For incompressible turbulence, the statistical interplay between temperature and velocity 
is decoupled, and thus, the temperature is regarded as a passive field ( Sreenivasan||1991 


Shraiman & Siggia 2000; Warhaft 2000). Nevertheless, for convective or compressible tur¬ 
bulence, because of the substantial impact to velocity through buoyancy or pressure, the 
temperature always behaves as an active field. Therefore, the properties of temperature 
statistics are the central issues. 


Belmonte & Libchaber (1996) experimentally studied the small-scale features of tem¬ 


perature in a thermal convection. They found that the temperature acts actively when 
the skewness product of temperature and its temporal derivative is positive, otherwise, it 
acts passively. The comparative experiments on the turbulent Rayleigh-Bernard convec¬ 
tion performed by Zhou & Xia (2008) showed that when the dynamical timescale is above 


(below) the Balgiano timescale, the temperature behaves as an active (passive) field. A 
fascinating feature of passive field is that the scaling exponent of structure function satu¬ 
rates for high order numbers, which is believed to be related to the so-called ”ramp-cliff” 
structures. Therefore, it is natural to ask whether a similar saturation appears for tem¬ 


perature. In fact, the findings in both experiments (Zhou & Xia 2002) and simulations 
(Celani et al. 2001; Zhou 2013) showed that the temperature is more intermittent than 
a passive field, and possesses a saturated exponent of 0.8, even smaller than the Burgers 
saturated exponent of 1.0 ( Mitra et al.| [2005). 

In terms of temperature in compressible turbulence, although it does not satisfy the 
standard definition of an active field, the nonlinear coupling to velocity makes it act 
actively. However, so far there are very few studies addressed this topic. In his theo¬ 


retical analysis, Canuto (1997) developed a model for handling compressible convection 


in the presence of large-scale flows, and obtained the dynamic equations for the mean 


and variance of temperature. Ni & Chen (2012) numerically investigated the statistics 


of temperature in one-dimensional compressible turbulence. They found that the tem¬ 
perature undergoes downscale cascade and follows the Kolmogorov picture. Moreover, 
the scaling exponent of temperature structure function is close to the Burgers scaling, 


indicating saturation at high order numbers. Recently, Donzis & Jagannathan (2013) 


carried out compressible turbulence simulations spanning the range of M t = 0.09 ~ 0.61. 
Their results were that: (1) the temperature spectrum defers to the fc -5 / 3 power law; 
(2) temperature fluctuations are less correlated to other thermodynamic variables, and 
the covariance between density and temperature contributes to the scaling of the mean 
and variance of pressure; and (3) the p.d.f. of temperature fluctuations is basically inde¬ 
pendent of turbulent Mach number. In detail, The p.d.f. tails for the positive component 
of temperature fluctuations are log-normal, while those for the negative component of 
temperature fluctuations retain a Mach-number dependence. 

Previous simulations in compressible turbulence ( Wang et a7||2011 2012a|b 2013a|b ) 
showed that there exists a strong connection between the shock topology and the forcing 
scheme. In particular, the flows driven by the solenoidal and compressive forces generate 
the small-scale shocklets and large-scale shock waves, respectively. This in turn, greatly 
influences the statistics of fields such as density, velocity and temperature. 

In this paper, we use the two groups data from the numerical simulations driven by the 
large-scale, random, solenoidal and compressive forcings to study the temperature in com¬ 
pressible turbulence, focusing on the effects of shock topology on the small-scale statistics 
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of temperature. The flows are computed on a 1024 3 grid by adopting a high-precision hy¬ 
brid method (Wang et al. J2010]), and the stationary values of the turbulent Mach number 
M t and Taylor microscale Reynolds number Re\ are (1.03, 255) for SFT and (0.62,164) 
for CFT, respectively. A systemic investigation on the fundamental statistics of tem¬ 
perature including the spectrum and field structure is actualized. Then, we report the 
statistics of dilatation, the application of isentropic assumption to the thermodynamic 
variables, and the angle statistics of temperature gradient on local flow structures. By 
employing a ” coarse-graining” approach to the temperature variance budget, we analyze 
the cascade of temperature, in particular, the crucial role of pressure-dilatation in the 
transport of temperature fluctuations. This paper is part of a series of investigations 


2015a 


on scalar transport in compressible turbulence. In three companion papers (Ni et al. 
Ni|20l5blct , we have examined carefully the statistical differences between active 


and passive scalars, the effects of Mach and Schmidt numbers on scalar mixing. We hope 
that this comprehensive study will advance our understanding of the small-scale statistics 
of temperature in compressible turbulence. 

The reminder of this paper is organized as follows. The governing equations and system 
parameters, as well as the details of simulation method, are described in Section 2. The 
basic statistics of the simulated flows is reported in Section 3. In the following two 
sections, we discuss the isentropic approximation of thermodynamic variables, and the 
statistical properties of temperature gradient. The analysis of the cascade of temperature 
is presented in Section 6 . In Section 7, we present the summary and conclusions. 


2. Governing Equations and System Parameters 

We simulate the statistically stationary compressible turbulence driven by the large- 
scale velocity forcing. Besides, a cooling function is added at large scales for removing 
accumulated internal energy at small scales. Similar to Ni et al. ( |2015a l, here we use the 
reference length L, density /?, velocity U and temperature Tq to normalize the compress¬ 
ible flow. Then we obtain the reference Mach number M = U/cq, where Cq = 1 /yRTg 
is the reference sound speed, R = C p — C v is the specific gas constant, and 7 = C p /C v 
is the specific heat ratio, with C p and C v representing the two specific heats at con¬ 
stant pressure and volume, respectively. By adding the reference dynamical viscosity 
p ,0 and thermal conductivity «o, we obtain another two basic parameters: the reference 
Reynolds number Re = UL/po and the reference Prandtl number Pr = poC p / kq- In our 
simulations, the values of 7 and Pr are set as 1.4 and 0.7, respectively. 

Based on the above procedure, the governing equations of the simulated flows, in 
dimensionless form, are written as 


d d , , 

p + rp\p u p = °> 


dt 

d 

dt 

d 

dt 


dxj 

d 1 

(pUi) + — [puiUj + pSij/'yM 2 ] = 


Re d Xj aij + pTi ' 


£+ + phM 2 ) Uj ] = - (n dT 


dxj 


a dx.j ' dx-j 


Re dxn 


P = pT, 


( 2 . 1 ) 

( 2 . 2 ) 

A -f- ppjUjj (2.3) 
(2.4) 


where a = PrRefiy—l)M 2 . The primary variables are density p, velocity it*, temperature 
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Table 1. Simulated parameters and resulting flow statistics. 


Flow 

Grid 

Re\ 

M t 

f^max^] 

T 

Lf L T f 

Ek 

Et 

e'/uj' 

SFT 

1024 3 

255 

1.03 

3.34 

1.11 

1.44 0.80 

2.13 

0.51 

0.34 

CFT 

1024 3 

162 

0.62 

3.18 

1.23 

1.56 1.12 

2.67 

0.53 

1.36 

S 3 

Ssm 

S3T 

u' T' 

p' 

u'c/u's 

<e> ( e *>/(e> 

<£c)/(e> 


(er) 

-2.2 

- 2.0 

0.18 

2.22 0.12 

0.29 

0.24 

0.54 88% 

14% 

-2% 

0.05 

-16.9 

-15.7 

0.76 

2.19 0.23 

0.52 

1.56 

0.79 32% 

68% 

0% 

0.16 


T and pressure p. Fj is the dimensionless large-scale velocity forcing 


Fj = X^i(ki)exp(*k/x) 


+ c.c.. 


(2.5) 


i=i 


Here i = y/— T, and Fj is the Fourier amplitude, which has only a solenoidal component 
perpendicular to k; for SFT ( Ni et al~[ 2013 2015a), but has another compressive compo¬ 
nent parallel to k; for CFT (Nij |2015 b). The viscous stress oy, and total energy per unit 
volume £ are defined by 


£ = 


f flu j 

K dxj 


P 


duj ■ 
dxi ■ 


1 


p( 


g pddij , 

( 2 . 6 ) 


(2.7) 


(7 — 1)7 M 2 2 

where 9 = duk/dxk is the velocity divergence or dilatation, a variable that quantifies the 
local rate of rarefaction (6 > 0) or compression (6 < 0). We now give the expressions 
of the dimensionless dynamical viscosity and thermal conductivity (Sutherland 1992), to 
complete the system 

1.4042T 15 . . 

tl,K= T +0.4042 ' (2 - 8) 

The large-scale velocity forcing presented in Equation (2.5) injects the same amount of 
energy into the two lowest spherical wavenumber shells. In particular, the energy injection 
in SFT is only perpendicular to the wavenumber vector, while that in CFT is both parallel 
and perpendicular to the wavenumber vector, and the ratio is 1 : 20. In addition, to keep 
temperature staying in a statistically stationary state, the velocity forcing and cooling 
function should satisfy the relation: (A) = ( pFjUj ), where (•) indicates ensemble average. 

The system is solved numerically by adopting a hybrid method in a cubic-box grid 
resolution of 1024 3 with periodic boundary conditions. This method applies a seventh- 
order weighted essentially non-oscillatory (WENO) scheme (Balsara & Shu 2000) for 
shock region and an eighth-order compact central finite difference (CCFD) scheme (Lele 
1992) for smooth region outside shock. A flux-based conservative formulation is imple¬ 


mented to optimize the handling at the interface between the shock and smooth regions. 
In addition, a new numerical hyperviscosity formulation is used to improve the numer¬ 


ical instability without compromising the accuracy of the hybrid method (Wang et al. 


2010). To obtain statistical averages of interested variables, a total of ten stationary flow 


fields were used. The values of large-eddy turnover time r are 1.11 in SFT and 1.23 in 
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CFT. Furthermore, in this paper, we always use the turbulent Mach number M t and 
Taylor microscale Reynolds number Re\, rather than the reference Mach number M and 
Reynolds number Re , to measure the compressible flows, which are defined as follows 


M t = 

(Vry 


Re\ = Re 


u'Mp) 
V3 (n)’ 


(2.9) 


where u' = ^{u\ + u 2 + u 3 ) is the r.m.s. velocity magnitude, and A is the Taylor mi¬ 
croscale, defined as A = u'/y/((dui/dxi) 2 + ( du^/dx-i) 2 + (du 3 /dx 3 ) 2 ). In our simula¬ 
tions, the stationary average values of ( M t , Re\) are (1.03, 255) for SFT and (0.62, 164) 
for CFT. 


3. Fundamental Statistics of Compressible Flows 

In this section, we first describe the simulated parameters of the compressible flows, 
then report the spectrum and field structure. The section is ended up with the discussion 
on the statistics of dilatation. 


3.1. Simulated Parameters, Spectra and Structure Fields 
Table [l] summarizes some overall statistics of the compressible flows. The resolution 
parameter k max r) are 3.34 in SFT and 3.18 in CFT, where p = [(p/(Rep)) 3 / < e/p >] 1//4 
is the Kolmogorov scale, and k max = N/2 = 512 is the maximum wavenumber in our 
simulations. It means that for both flows, the fine-scale structures in smooth regions are 
well resolved by the hybrid method. Here we point out that although the thickness of 
shock is comparable to the grid length and is not directly resolved by the WENO scheme, 
the total amount of dissipation across shock is independent of numerical viscosity. 

The integral scales for velocity and temperature are computed by 


3tt /' E(k) 7T f E T (k) 

Lf = 2^ 2 ~JT dk ’ LT f = T^ ~l^ dk ' 


(3.1) 


where E{k) and ET{k) are the spectra of kinetic energy and temperature per unit mass, 
respectively 


J E(k)dk = ^u' 2 , J E T (k)dk = 


(3.2) 


o o 

The ratio of the r.m.s. magnitudes of dilatation to vorticity, d'/ui', are 0.34 in SFT and 
1.36 in CFT, showing that at small scales, the effect of compressibility in is overwhelming 
in the compressive forced flow. Here 9' = \J (9 2 ) and u/ = \J (w 2 + uj 2 + cc|) are the r.m.s. 
magnitudes of dilatation and vorticity, respectively. 

The compressible character can be further demonstrated by ensemble averages of the 
skewness of velocity derivative and the mixed skewness of velocity-temperature deriva¬ 
tives 

[<(&)’ +<fe) 3 + (fe) S }]/ 3 (33) 


S, = 


[<(fe ) 2 + <! 5> 2 + ( 


$3 m — 


f / dui r dT \2 

[A fa*i > 


du 2 

dx 2 


( dT 

' 6X2 


dx 3 
du 3 


) 2 >/ 3 ]' 


\2 i du 3 / 
> ^ dx 3 V 


,H) 2 )]/3 


[<(fe ) 2 + (fe ) 2 + (It) 2 )/ 3 ] /_ [((Ih ) 2 + (H ) 2 + (H) 2 )/ 3 ] 


(3.4) 


... dT ^2 j ( dT ' 

'dx 2 ^ 1 ^ dx 3 > 

In SFT S 3 and S m3 are —2.2 and —2.0, respectively; however, in CFT their magnitudes 
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increase to S 3 = —16.9 and S 3 rn = —15.7. This indicates that the presence of large-scale 
shock waves leads the velocity field of CFT to be anisotropic. By contrast, for ensemble 
average of the skewness of temperature derivative 


S3T 


[«g ) 3 + (l £) 3 + (IS) 3 )]A 

trd 0.76 for CFT, implying that 


(3.5) 


The values are 0.18 for SFT and 0.76 for CFT, implying that the temperature field 
is approximately isotropic in compressible turbulence, espcially in the solenoidal forced 
flow. In addition, the r.m.s. magnitude of temperature fluctuations is 0.12 in SFT and 
0.23 in CFT, while that of density fluctuations reaches as high as 0.29 in SFT and 0.52 
in CFT. This reveals that compared to temperature, the compressibility makes density 
vary more intense. 


The application of Helmholtz decomposition (Samtaney et al. 2001) to velocity field 
gives that 

u.i — Uis + ( 3 . 6 ) 

where Ui S is the solenoidal component satisfying dui S /dxi = 0 , and Ui C is the com¬ 
pressive component satisfying £ijkduk c /dxj = 0, with representing the Levi-Civita 
symbol. The ratio for the r.m.s. magnitudes of the two components, u' c /u' sl is 0.24 in 
SFT and 1.56 in CFT. Furthermore, Andreopoul os et aZ.| (2000) showed that the vis¬ 
cous dissipation rate e = (JijSij/Re can be divided into three parts: the solenoidal 
part e s = (fj,/Re)uJiUii, the dilatation part e c = (4/3)(/j,/i?e)0 2 , and the residual part 
e m = (2 n/Re)[(dui/dxj)(duj/dxi) — 9 2 ]. In our simulations, the percentages of the first 
two parts are ( 88 %, 14%) for SFT and (32%, 68 %) for CFT, meaning that in SFT most 
kinetic energy is dissipated through vortices stretching, while in CFT the dissipation of 
kinetic energy is dominated by rarefaction and compression.The temperature dissipation 
rate is defined as 

e T = n(dT/d Xj ) 2 . (3.7) 

Table |T] shows that the ensemble-average value, (er), is 0.05 in SFT and 0.16 in CFT, in¬ 
dicating that compared to small-scale shocklets, there are more temperature fluctuations 
depleted by large-scale shock waves. 

Figure [T] presents the compensated kinetic energy spectra from the two simulated 
flows. Here an operational definition of the inertial range is identified for the overall 
kinetic energy spectrum, namely, Ck = E(k)k 5 / 3 /(e) 2 / 3 , where Ck is the Kolmogorov 
constant. In SFT, Ck is found to be 2.17 for the range of 0.03 < kr] < 0.1, a bit higher 


than the typical values observed in incompressible turbulence (Wang et al. 1996). By 


contrast, for the same range is CFT, Ck has a much smaller value of 1.05. In the inset 
we plot the solenoidal and compressive components of the kinetic energy spectra. For the 
solenoidal components, the appearance of spectral bumps at high wavenumber reveals 
the nonlocal feature for the transfer of the solenoidal component of velocity fluctuations 
in the crossover region between inertial and dissipative ranges. Moreover, it shows that 
the compressive component in SFT follows the fc -5 / 3 power law, while that in CFT defers 
to the k ~ 2 power law over than one decade wavenumber range. As one known, previous 


studies on Burgers turbulence (Bee & Khanin 2007) showed that it is the existence of 
large-scale shock wave gives rise to the k ~ 2 scaling. 

In Figure [2] we plot the compensated spectra of temperature. Similar to that in com¬ 
pressible flow, the temperature spectrum in SFT displays the k ~ 5 / 3 power law, and its 
scaling constant can be operational defined as Ct = ET(k)k 5 ^ 3 (e ) 1 ^ 3 / (er) ■ It shows that 
in the range of 0.06 ^ ^ 0 . 12 , the value of Ct is about 0.08, much smaller than the 


















7 


Effects of shock topology on temperature field in compressible turbulence 



Figure 1. Compensated spectrum of kinetic energy, where the dashed and solid lines are for 
SFT and CFT, respectively. Inset: solenoidal (solid) and compressive (open) components, where 
the symbols of squares and circles are for SFT and CFT, respectively. 



Figure 2. Compensated spectrum of temperature, where the dashed and solid lines are for SFT 
and CFT, respectively. The dotted line is for the compensated spectrum of density-weighted 
temperature in CFT. 


typical values observed from passive scalar. By contrast, the spectrum of temperature in 
CFT has the k~ 2 power law, indicating the dominant motion of large-scale shock wave. 
In the inset we present the spectrum of the density-weighted temperature, p 1 ^T , from 
CFT. The result is that the spectrum still follows the k~ 2 power law, which is different 
from the k~ 5 / 3 power law for the spectrum of the density-weighted velocity, p 1 ^ 3 u, re¬ 
ported by Kritsuk et al. (2007). It implies that in compressible turbulence the density 
and temperature do not have the direct statistical coupling, they only connect with each 
other through velocity. 

We now take attention on the field structures. Figure [3] provides the two-dimensional 
(2D) contours of temperature (top left), temperature dissipation rate (top right), vortic- 
ity magnitude (bottom left) and dilatation (bottom right) in SFT. Here the temperature 
field has the ”ramp-cliff” structures, to some extent, is like the scalar field in incompress¬ 
ible flows (Shraiman & Siggia 1994). In particular, the small-scale cliffs with high gradi- 
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Figure 3. Two-dimensional contours of temperature (top left), temperature dissipation rate 
(top right), vorticity magnitude (bottom left), and dilatation (bottom right) in SFT. 


ents of fluctuations divide the large-scale ramps with low gradients of fluctuations. The 
temperature dissipation field consists of the smooth, low-dissipation sea and the sharp, 
high-dissipation discontinuities under random distribution. For the dilatation field, the 
sharp discontinuities are exactly the small-scale shocklets. Moreover, the intermittency of 
the vorticity field is reflected by the inhomogeneous distribution of the small-scale strong 
vortices. The same four 2D contours of CFT are depicted in Figure [4j There appear 
large-scale discontinuities in the temperature, temperature dissipation and dilatation 
fields, causing by the large-scale shock waves. In front of shock waves, the temperature 
fluctuations are small, however, they increase quickly once across the shock waves. It 
further shows that for the vorticity field, the small-scale strong vortices preferentially 
concentrate in the high temperature regions, and thus, displays intensive intermittency. 


3.2. Probability Distribution Functions 

In Figure [5] we plot the p.d.f.s of the normalized vorticity magnitude. As expected, they 
exhibit well-defined exponential tails, and the one in CFT is very long, showing strong 
vorticity fluctuations in the condition that the compressive mode of velocity is stimulated. 
Here we also present the incompressible p.d.f. provided by Moisy & Jimenez (2004). At 
large amplitudes, the p.d.f. collapses between those from SFT and CFT. This reveals 
that the declaration in Wang et al. (2012b) that the intense vorticity was suppressed in 
compressible turbulence works only for the solenoidal forced flow. 

The log-log plot of the p.d.f.s against the normalized negative dilatations is shown in 
Figure [6] In strong compression region, the p.d.f.s have the power-law tails, which are 


qualitatively similar to the p.d.f. of velocity derivative in Burgers turbulence (Bee & 
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Figure 4. Two-dimensional contours of temperature (top left), temperature dissipation rate 
(top right), vorticity magnitude (bottom left), and dilatation (bottom right) in CFT. 



Figure 5. The p.d.f. of normalized vorticity magnitude, where the dashed and solid lines are 


for SFT and CFT, respectively, and the circles are from incompressible turbulence (Moisy & 
Jimenez|2004 1 


Khanin|2007). According to the stochastic theory of Burgers equation, it is the preshock 


leads to the large negative velocity derivative, and then the power-law tail (E et al. 1997 E 
fe Eijnden|1999 2000 Bee et a/.|2000[ Bec|2001 ). In our simulations the exponents of the 


power-law tails are around —2.5 in SFT and —3.5 in CFT. We notice that the exponent 
value in the compressive forced flow is the same to that in the Burgers turbulence. Indeed, 
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Figure 6 . Log-log plot of the p.d.f. of normalized dilatation in compression region, where 
the dashed and solid lines are for SFT and CFT, respectively. Inset: the p.d.f.s of normalized 
dilatation. 


previous study (Wang et al. 2012a) showed that the power-law region in compressible 
turbulence is mainly caused by preshocks and weak shocklets rather than strong shock 
waves. In the inset we plot the p.d.f.s throughout the entire dilatation range. Obviously, 
they exhibit strong skewness towards negative side. 

To understand the mechanism of the power-law behavior, we write the Liouville equa¬ 
tion for the p.d.f. of dilatation, Q(9), as follows 


d_ 

dt 


Q 


d_ 

d6 


(P + R- D + 9 2 )Q 


9Q. 


(3.8) 


Here P(6), R(9 ) and D{9) are the ensemble averages of pressure, anisotropic straining 
and viscous dissipation conditioned on dilatation, respectively 


p <*>=C4 ( ^ )|9> ' 

(3.9) 


(3.10) 


(3.11) 


In Figure [7] we first plot the conditional averages of straining and its compressive com¬ 
ponent, S(9/9') = {duj/dxidui/dxj/9' 2 \6/0') and S c (9/9') = ( du’//dx i du'//dxj/9' 2 \9/9'), 
as functions of the normalized dilatation 9/9'. It shows that for in the compression re¬ 
gions for both flows, the compressive component dominates the overall straining. Namely, 
the contribution from the solenoidal component is negligible. We then plot P{9/9') and 
D(9/9') against 9/9'. Similar to Wang et al. (2012b), in SFT P{9/9') and D{9/9') are 
positive and can be well described by a quadratic-parabolic formulation of C\ (9/9') 2 + 
C2{9/9'). The values of ci and C 2 are the same to those in Wang et al. (2012b). They 
are c\ p = 1.2, C 2 P = 5.3, C\d = 0.4 and C 2 d = 0.5, where the subscripts p and d are for 
pressure and viscous dissipation, respectively. By contrast, in CFT the values of C\d and 
C 2 d are 0.2 and 0.4, respectively. However, P{9/9') is negative, and its fitting coefficients 
are Ci p = —0.1 and C 2 P = 0.3. 
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Figure 7. S(9/9') (circles), S c (6/6') (deltas), P(9/9') (diamonds) and D(9/9') (circles) as func¬ 
tions of normalized dilatation. The lines are for the fitting by a quadratic-parabolic formulation 
at various values of coefficients, (a) SFT, (b) CFT. 



Figure 8. The p.d.f.s of normalized temperature increment, where asr is the standard 
deviation of temperature increment, (a) SFT, (b) CFT. 


For large negative dilatation, the stationary solution of Equation (3.7) 

Q(0) oc Q~fi 


(3.12) 


where q = 2+l/(l+ci p —cid). Immediately, we obtain that q = 2.56 for SFT and q = 3.43 
for CFT. which are very close to the exponent values of 2.5 and 3.5 shown in Figure [6j 
It is worth to point out that in SFT the contributions to the power-law exponent from 
the pressure and viscous dissipation are converse, while in CFT they become the same 
because of the reversal of the contribution from pressure. 

In Figure[ 8 ]we plot the two-point p.d.f.s of the normalized temperature increment at the 
normalized separation distances of r/rj = 1, 4, 16, 64 and 256. Here ST = T(x + r) — T(x) 
is the temperature increment. It shows that the p.d.f. is basically symmetric for each 
scale, and gets broader as scale increases. In particular, in SFT the p.d.f. at r = 256?; is 
approximately Gaussian, while that in CFT still keeps super-Gaussian. This implies that 
the presence of large-scale shock wave makes the temperature field be intermittent, even 
though it lies on the integral scale. At large amplitudes, the shape of p.d.f. is concave, 


similar to that of longitudinal velocity increment shown in Ni et al. (2015a). 


The rescaled p.d.f., Q(ST/T l )/r , for the inertial range of r/rj = 16, 32 and 64 are 
shown in Figure [9] It is found that in CFT the p.d.f.s collapse to the same distribution. 
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Figure 9. The rescaled p.d.f.s of normalized temperature increment, (a) SFT, (b) CFT. 



Figure 10. T/(T) versus (p/ip)) 1 where 2197 points are used in the plot. The dashed 
line is for isentropic relation, (a) SFT, (b) CFT. 


According to multifractal theory (Benzi et al. 2008), the scaling exponents for the sta¬ 
tistical moments of temperature increment in the inertial range of CFT will saturate at 
high order numbers. In other words, (|AT| P ) oc r z °° for large p , and zoo = 1. 


4. Isentropic Approximation on Thermodynamic Fields 

In compressible turbulence it is often valuable to reduce the number of thermodynamic 
variables, for simplifying theoretical analysis and developing engineering model. A wide 
used approach, though not strictly exact because of the irreversible dissipative nature 


of turbulence, is to assume that thermodynamic processes occur isentropically (Chan¬ 


drasekhar 1951 Erlebacher et al. 1990). According to the state equation, the conservation 


of energy at constant entropy indicates that the instantaneous density, temperature and 
pressure are related as follows 


t/(t) = ( p /(p)y 1 = (p/<p» 


V 1—1/7 


(4.1) 


In Figure 10 we plot T/(T) against ( p/(p)) 1 1/7 . A total of 2197 points are presented, 


and the dashed line stands for the isentropic relation. In SFT most of the points collapse 
onto the isentropic line. There are only about 1.6% points failing to satisfy the isentropic 
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Figure 11. Conditional averages of ratio of T/(T) to ( p/{p}) 1 1 '~ r , where the dashed and solid 

lines are for SFT and CFT, respectively. 


relation within a tolerance of Ctoi = 10%, which is defined by 

Ctoi = |1 - {p/{jp)) l ~ lh /{T/{T))\. 


(4.2) 


Nevertheless, within the same tolerance there are about 18.9% points failing to satisfy 
the isentropic relation in CFT. The large-scale shock waves make the fluctuations of ther¬ 
modynamic variables intense and destroy the insentropity, especially at large amplitudes. 
We then define the correlation coefficient between T/(T) and (p/ip)) 1-1 ^ 1 as follows 


Ct p = 


(T/(T) — l) ((p/(p)) 1-1 / 7 — l) 


(T/(T)-i) 2 \ 1/2 (((p/( P ))i^/7 


W> 


1 / 2 ' 


(4.3) 


It shows that the values of Ct p for both SFT and CFT are 0.96. 

The average of the ratio of T/(T) to (p/(p)) 1-1 ^ 7 conditioned on dilatation, as a 


function of the normalized dilatation, is shown in Figure 11 In the compression region 
of —10 < 9/9' < —2, the average value in SFT is near unity, and slightly decreases as 
compression decreases. This implies that the isentropic approximation of thermodynamic 
variables is valid in the small-scale shocklet regions. By contrast, the average value in 
the same compression region of CFT is a bit higher, and is close to 1.04 in the range 
of —8 ^ 9/9' ^ —2.4. In a word, for the compression region of compressible turbulence, 
one can use the isentropic approximation to facilitate the description of thermodynamic 
variables. 

As one known, the thermodynamic fluctuations are related through Equation (4.1) on 
an instantaneous basis. We now introduce a quantity of 

P/(P) 


_ P/(P) 

tp = 


(. P/(P)Y ( T/(T )) ; 


(4.4) 


to measure the exactitude of isentropic approximation. Obviously, if the case is isentropic, 
the distribution of ip' = <p — (ip) will follow the Dirac delta function, Q(p') = S(p'). In 


Figure 12 we plot Q{<p') from SFT and CFT. For comparison, those at M t = 0.1, 0.3 and 


0.6 in Donzis & Jagannathan (2013) are also presented. It shows that in the solenoidal 
forced flow with M t = 0.1, Q(<p') is very close to 6(<p'), and deviates from that as M t 
increases. At M t « 0.6, Q(p') of CFT is much wider than that of SFT. Furthermore, 
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Figure 12. Q(p') in SFT (squares) and CFT (circles).The d otted, dashed an d solid lines are 
for those at M t = 0.1, 0.3 and 0.6 flows presented in Donzis & Jagannathan (2013). 


the value of the standard derivation, ((tp' / (ip)) 2 ) 1 / 2 , which quantifies the departure from 
isentropic assumption, is 0.19 in SFT and 0.27 in CFT. Here we point out that the results 
from our simulations do not follow the 0.1 M 2 scaling proposed by 
((2013). 

The above analysis shows that globally, the thermodynamic variables in compressible 
turbulence are not exactly isentropic, because of Q{<p') deviates from 6(<p'). Nevertheless, 
by introducing an exponent parameter 7 C , we can still use the following expression to 
connect pressure and temperature 

p/(p) = (T/(T))^. (4.5) 

It gives that 

tp = (T/iT))^-^. 

We then use the constraint of probability, Q(ip)dip = Q(T)dT , to evaluate to what 
7 C can capture the fluctuations of thermodynamic variables. It yields that 

Q{p) = \ (T)(T/(T)) 1+ ^-^Q((T)p^~^). (4.7) 

I 7c — 1 7—1 I 

By assuming that the temperature field follows Gaussian distribution, we obtain the 
expression of Q(ip) 


(4.6) 

degree 


Donzis & Jagannathan 


7 — 1 7c —1 


Q(<p) = 


p 


V2tt(T* 2 )i/ 2 


■ exp 


2 (T* 2 ) 


/ -J. _\ 5 

[ip-t- 1 Tc-1 — II 


(4.8) 


where T* = T / (T) is the normalized temperature fluctuation, with T representing the 
temperature fluctuation. 


In Figure 13 we compare Q(p) between theoretical model and numerical simulation. 

can be obtained through the best fitting of simulation data, 
in SFT and CFT are around 1.20 and 1.22, respectively, 


For a certain condition, 7 , 

In our cases, the values of 7 , 
smaller than the isentropic value of 1.4. It displays that Q(<p) from SFT is basically well 
described by the model, except at the large positive amplitudes. By contrast, Q(<p) from 
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Figure 13. Q(ip) in SFT (squares) and CFT (circles). The dashed and solid lines represent 
the model results for SFT and CFT, respectively. 



Figure 14. The p.d.f.s of cosine of angle between temperature gradient and vorticity, where 
the dashed and solid lines are for SFT and CFT, respectively. 


CFT goes far away from the model, which indicates that the Gaussian assumption of 
temperature field in CFT works badly. 


5. Statistical Properties of Temperature Gradient 

Contrasting to that the properties of temperature in turbulence have received consid¬ 
erable attention, there are few studies addressed the properties of temperature gradient. 
In this section we study the statistics of temperature gradient on local flow structures. 
In Figure |l4| we plot the p.d.f. of the cosine of the angle between temperature gradient 
g t = (dT/dxi)/(dT/dxi) 2 + ( dT/dx 2 ) 2 + ( dT/dx 3 ) 2 and vorticity. It shows that in 
highly compressible turbulence, there is a strong tendency for the temperature gradi¬ 
ent being orthogonal with the vorticity, especially for CFT. This feature is similar to 
the observation from passive scalar transport in turbulence (Ni et al. 2015a). Thus, we 
speculate that in some respects, temperature in compressible turbulence may behave like 
passive scalar when the degree of compressibility increases. 
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Figure 15. The p.d.f.s of cosines of angles between temperature gradient and vortex stretching 
vector as well as its components, where the solid, dashed and dotted lines are for the vortex 
stretching vector and its solenoidal and compressive components, respectively, (a) SFT, (b) 
CFT. 


As a further investigation, we now consider the statistical correlation between the tem¬ 
perature gradient and the anisotropic strain rate tensor S*j — Sij — Skk$ij/3 (Erlebacher 


| et a/.|1993 Pirozzoli & Siggia 2004), which is written through the enstrophy equation as 
follows 


d 

Ft ' 

+ UJi 


d ' 

^ dxj / 




■uiu j s%--ew 2 


£ ijk 1 dp dp 
7 M 2 p 1 


dxj dxk 


■ Uli 


£ijk 


d /l d<Jmk \ 
dx.j V p dx m / 


Re 

j -r duj/dxj)/2, and Sfj and Sf 
respectively. Further, the vortex stretching vector W, = 


(5.1) 


Here Sij = ( dui/dx 
parts of S* 

as Wf = oJjSfj and Wf = uijSf*. In Figure [l5| we plot the p.d.f.s of the cosines of the 


are the solenoidal and compressive 
is decomposed 


angles between temperature gradient and vortex stretching vector and its components. 
Globally, in SFT there are small positive alignments between the temperature gradient 
and the vortex stretching vector and its components. In detail, the maximum p.d.f. for the 
solenoidal component appears in the case where it is perpendicular to the temperature 
gradient, while that for the compressive component appears in the case where it aligns 
with the temperature gradient. By contrast, in CFT the p.d.f.s for the vortex stretch¬ 
ing vector and its components are basically symmetric. It shows that the temperature 
gradient is preferentially perpendicular to the solenoidal component, and preferentially 
aligns with the compressive component. These results reveal that in compressible tur¬ 
bulence the increase in the degree of compressibility suppresses the anisotropy of strain 
rate tensor. 

Let us denote the three eigenvectors of the anisotropic strain rate tensor as Ai, A 2 
and A3. The corresponding eigenvalues, arranged in ascending order, are A*, A 2 and A 3 , 
which satisfy the following condition 


At + A 2 + A 3 = 0. (5.2) 

Figure [16] presents the alignment statistics between the temperature gradient and the 
strain rate eigenvectors. The results are that: (1) there is a strong tendency for the tem¬ 
perature gradient to align with the first eigenvector corresponding to the most negative 
eigenvalue; ( 2 ) there is also a clear tendency for the temperature gradient to be per¬ 
pendicular to the second eigenvector; and (3) the tendency for the temperature gradient 

















Effects of shock topology on temperature field, in compressible turbulence 


17 




Figure 16. The p.d.f.s of cosines of angles between temperature gradient and the eigenvector Aj, 
where i = 1,2 and 3. The lines with open and solid symbols are for SFT and CFT, respectively. 
Inset: the p.d.f.s for the solenoidal (dashed line) and compressive (dotted line) components. 


to be perpendicular to the third eigenvector, the one with the most positive eigenvalue, 
is also noticeable. The insets of Figure 16 show that for each flow the solenoidal com¬ 
ponent dominates the alignment statistics, and the contribution from the compressive 
component mainly occurs at small angles. 

The conditional averages of the squares of cos(g, A*) and are plotted in Figure [l7| As 
the compression increases, the conditional average for the first eigenvector approaches 
unity, while those for the second and third eigenvectors approaches zero. This indicates 
that in strong compression region, the temperature gradient always aligns with the most 
negative eigenvalue of the strain rate tensor. In fact, in the vicinity of a shock, the tem¬ 
perature gradient and the first strain rate eigenvector are both perpendicular to the 
shock front. Furthermore, in rarefaction region, there is a clear tendency for the tem¬ 
perature gradient to align with the third eigenvector corresponding to the most positive 
eigenvalue. 


Besides the angle statistics of temperature gradient, we now explore the effects on 
temperature gradient from local flow structures. To facilitate the description of local flow 
structures in compressible turbulence, we introduce the first, second and third invariants 
of the anisotropic velocity gradient tensor A*j = A.^ — dSij/3 as follows 


Wang et al. 
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Figure 17. Conditional averages of squares of cosines of angles between temperature gradi¬ 
ent and the eigenvector Ai, where the solid, dashed and dotted lines are for * = 1, 2 and 3, 
respectively. The open and solid symbols are for SFT and CFT, respectively. 



Figure 18. Isocontour lines of logarithm of joint p.d.f. of Q* and R*. Five contour lines at 0, 
— 1, —2, —3 and —4 are shown, (a) SFT, (b) CFT. 


2012 a) 


p* = -(& + & + &) = °’ 

q* = m+ ess+ =Q-\ p2 

R* = =R- 1 -pq+^p\ 


(5.3) 

(5.4) 

(5.5) 


where £* = & — 0/ 3 are the three eigenvalues of A*-, and & are the three eigenvalues of 


A,,j = duj/dxi . The details in P, Q and R can be found in Chong et al. (1990). 


In Figure [18] we plot the contour lines of the joint p.d.f. for the second and third 
invariants ( Q*,R *), which are scaled by the second invariant of the rotation rate tensor 
Qn = flijClij/2. Similar to the results from Pirozzoli fe Siggia| (2004), the joint p.d.f.s 
display the teardrop shape in both SFT and CFT. Compared with weakly and moderately 
compressible turbulence (Pir ozzoli fc Siggia| 2004),in our simulations the tails of the joint 
p.d.f.s in the fourth quadrant are pronounced longer. It also shows that the joint p.d.f. 
from CFT has protrusive structures in the third quadrant. Figure[l9|presents the contour 
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Figure 19. Isocontour lines of temperature gradient magnitude conditioned on Q* and R*. 
Nine lines at 0, 0.25, ..., 2 are shown, (a) SFT, (b) CFT. 




Figure 20. Subgrid-scale temperature flux (circles), pressure-dilatation (squares), and viscous 
dissipation (deltas), as functions of normalized scale l/ij. (a) SFT, (b) CFT. 


lines of the magnitude averages of temperature gradient conditioned on Q* and R*. There 
also appears the teardrop in the conditional average, which shifts towards the positive 
part of Q* and the negative part of R*. This reveals that the temperature gradient 
is relative large in the region where the second and third invariants of the anisotropic 
velocity gradient tensor are respectively positive and negative. 


6. Cascade of Temperature Field 

In incompressible flow, the cascade of temperature involves the generation of temper¬ 
ature fluctuations at large scales, the stretching, contracting and folding of temperature 
by velocity, producing progressively smaller and smaller scales, until the ultimate dissi¬ 
pation of temperature fluctuations at the smallest scale. In compressible turbulence, the 
fact of the nonlinear interplay between velocity and temperature as well as that between 
solenoidal and compressive modes of velocity greatly complicate the cascade of tempera¬ 
ture. In this section, we carry out investigation on this topic, especially in analyzing the 
important role of pressure-dilatation in temperature cascade. 

A ”coarse-graining” approach (Al uie|2011 Al uie et al.|2012 ) is employed to study the 
transport of temperature fluctuations at different scales. According to the definition of a 
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Figure 21. Normalized solenoidal (solid symbols) and compressive (open symbols) components 
of the subgrid-scale fluxes of temperature, as functions of normalized scale l/r), where the symbols 
of squares and circles are for SFT and CFT, respectively. 


classically filtered field <q(x) 


az(x) 


d 3 rGi(r)a(x + r), 


( 6 . 1 ) 


the density-weighted filtered field is obtained 


a*(x) 


PQz(x) 
dz( x ) ’ 


( 6 . 2 ) 


Here G;(r) = G(r/l)/l 3 is a kernel, and G(r) is a window function. By the large-scale con¬ 
tinuity and temperature equations, it is straightforward to derive the governing equation 
for temperature variance at large scales as follows (Ni et aZ.||2015al 


d ( !-^r 2 
dt\2 P ' Te ‘ 


V • J I = —n; — <bz — A; — Di + £' 


cool 


(6.3) 


In Figure |20| we plot the ensemble averages of SGS temperature flux, pressure-dilatation, 
and viscous dissipation, as functions of the normalized scale l/r]. In both SFT and CFT, 
the viscous dissipation primarily occurs at small scales, and declines quickly as scale 
increases. Throughout scale ranges, the SGS temperature flux is positive, indicating 
that the temperature fluctuations are always transported from large to small scales. 
The appearance of plateau in the SGS temperature flux confirms the conservation in 
temperature cascade. In terms of the magnitude of pressure-dilatation, in SFT it first 
increases and undergoes a flat region in the range of 8 < l/r] < 24, then decreases 
and reaches zero at large scales. By contrast, in CFT the flat region shifts to a larger 
scale range of 16 ^ l/r] < 80, and when scale increases it decreases and approaches 
a finite positive value. These observations reveal that in the transport of temperature 
fluctuations, the pressure-dilatation mainly takes activities at moderate and large scales. 

The solenoidal and compressive components of the SGS temperature flux normalized 
by the ensemble average of temperature dissipation rate are depicted in Figure [21] Ob¬ 
viously, the two components are positive throughout scale ranges. In SFT the solenoidal 
component is smaller than the compressive component in the range of l/r] ^ 96, and 
becomes comparable at larger scales. By contrast, in CFT the solenoidal component is 
always smaller than the compressive component, revealing the overwhelming effect of 
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Figure 22. The p.d.f.s of cosine of angle between Tj(u,T) and dT/dxj (solid lines), and its 
solenoidal (dotted lines) and compressive (dashed lines) components, (a) SFT, (b) CFT. 


compressibility. To explore the angle statistics of temperature gradient in the cascade 
process, we rewrite the expression of the SGS temperature flux as follows 


II; = —Tj (u, T)dT/dxj , 

where the SGS temperature-velocity coupling fj is defined as 

t 3 {u,T) 

Then, we obtain the angle between fj and dT/dxj 


cos p = 


jdT/da 


i\\dT/da 


(6.4) 


(6.5) 


( 6 . 6 ) 


In Figure 22 we plot the p.d.f.s of the cosine of the angle (3 and its components (3 S and 
fi c . For both SFT and CFT, the p.d.f.s of cos j3 are asymmetric and peak at j3 = 0. This 
indicates that the transfer of temperature flux preferentially occurs in the orientation 
where fj is anti-aligns with dT/dxj. The p.d.f.s of cos f3 s are also asymmetric, but peak 
in the range of 0 < fi s < 7t/2. On the contrary, the p.d.f.s of cos /3 C are basically symmetric 
and peak at fi c = 0 and n. 

The pressure-dilatation PD(l ) = pfii plays an important role in the conversion between 
kinetic and internal energy, namely, if 9i < 0, the energy is transported from the large- 
scale kinetic to internal energy; if 9i > 0, the process reverses. In Figure 23 we plot the 
cospectrum of pressure-dilatation, which is defined by 


E pd (k) = Y, -P(k)0(-k). 

fc-l/2<|k|<fc+l/2 


(6.7) 


It shows that the slope value of the cospectrum is —1.30 for SFT and —2.38 for CFT, 
exhibiting that in our simulations E pd (k) decay at rates faster than fc -1 . This is in 
agreement with the following criterion proposed by Aluie (2011) 


\E pd (k)\ < Cu'p'(kL)~ 


w > 1, 


( 6 . 8 ) 


where C is a nondimensional constant and L is an integral scale. This criterion deduces 
that the pressure-dilatation would converge and become independent of l at small enough 
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Figure 23. Cospectrum of pressure-dilatation, where the dashed and solid lines are for SFT 

and CFT, respectively. 



Figure 24. Pressure-dilatation (circles) and its positive (deltas) and negative (squares) 
components, where the open and solid symbols are for SFT and CFT, respectively. 


scales 


limPDCO = mrio £ E pd (k) = e 

0 <k<K 


(6.9) 


It provides the picture that the pressure-dilatation mainly exchanges the kinetic and 
internal energy over moderately large scales of limited extent. At smaller scales the 
kinetic and internal energy budgets statistically decouple, giving rise to an inertial range 
in which the temperature undergoes a conservative cascade. Furthermore, in CFT the 
faster decaying of E p d(k ) undoubtedly connects with the larger compressive component 
of the SGS temperature flux. 

The pressure-dilatation and its positive and negative components, as functions of the 
normalized scale l/r], are depicted in Figure 24 Throughout scale ranges, and ($j~) 


themselves are substantial, however, when adding together they almost cancel each other 
and make the outcome small, especially at small scales. In other words, the high values 
of pressure-dilatation generated in the vicinity of shock structures will basically vanish 
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Figure 25. The p.d.f.s of the pressure-dilatation for large scales (a) and the residual part from 
small scales (b) at l = 32 g, where the dashed and solid lines are for SFT and CFT, respectively. 
The dotted lines are for Gaussian. 



Figure 26. Two-dimensional contours of the pressure-dilatation for large scales (a) and the 
residual part from small scales (b) at l = 32r/, in SFT. 


after taking global averages. Further, we point out that the picture of the negligible 
contribution of pressure-dilatation at small scales does not contradict to the motions of 
rarefaction and compression appearing at all scales, which is the fundamental property 
of compressible turbulence. In CFT (<f>^) is slightly larger than that of (<!>+}, while in 
SFT the tw o co mponents mostly overlay throughout scale ranges. 

In Figure 


25 


we plot the p.d.f.s of the pressure-dilatation for large scales pfii, and 
the residual part from small scales p9 — pfii- The results are that the p.d.f. of pfii is 
sub-Gaussian and has small positive skewness. Compared to those in CFT, the p.d.f. 
tails in SFT are broader and thus are less intermittent. The super-Gaussian p.d.f.s of the 
residual parts display heavy tails, implying spatially rare but intense two-way exchange 
between kinetic and internal energy. Moreover, these p.d.f.s are strongly positive skewed, 
indicating that the conversion from kinetic to internal energy through compression is 
more intense than the inverse process through rarefaction. In Figures |26| and [27} we plot 
the 2D contours of the pressure-dilatation for large and small scales. The filter width 
is l = 32 rj. Obviously, there exists appreciable differences in the behavior of pressure- 
dilatation between the solenoidal and compressive forced flows. In SFT the pressure- 
dilatation for large scales and the residual part from small scales distribute randomly, 
and the discontinuities shown in Figure [26^ b) are the small-scale shocklets. By contrast, 











24 


Qionglin Ni 



Figure 27. Two-dimensional contours of the pressure-dilatation for large scales (a) and the 
residual part from small scales (b) at l = 32 77 , in CFT. 



Figure 28. Fluctuation component of pressure-dilatation. The open and solid circles are for 
SFT and CFT, respectively, and the dashed and solid lines represent the results computed from 
the Zeman model. 


in CFT the two parts pressure-dilatation concentrate in the rarefaction and compression 
regions, and the discontinuities shown in Figure |27[ b) are the large-scale shock waves. 
In Figure [28] we depict the fluctuation component of pressure-dilatation, as a function of 
the normalized scale Z/ 77 . For comparison, we also present the results computed from the 
Zeman model ( Zeman[ 1991). It shows that (p'9) from SFT basically collapses onto the line 
representing the model, except a few derivations at small and large scales. Nevertheless, 
because of the strong degree of compressibility, (p'9) from CFT deviates far away from 
the model. 


7. Summary and Conclusions 

In this paper, we performed a systematic investigation on the effects of shock topology 
on the statistics of temperature in isotropic compressible turbulence. The simulations 
were solved numerically using a hybrid method of a seventh-order WENO scheme for 
shock region and an eighth-order CCFD scheme for smooth region outside shock. The 
small-scale shocklets and large-scale shock waves appeared in the compressible turbulence 














25 


Effects of shock topology on temperature field in compressible turbulence 

driven by the large-scale, solenoidal and compressive forcings, respectively, where the 
stationary values of the turbulent Mach number and Taylor microscale Reynolds number 
were (M t , Re\) = (1.03, 255) for SFT and (M t , Ref) = (0.62, 164) for CFT. A variety of 
issues including the spectrum, field structure, probability distribution function, isentropic 
assumption and cascade were reported. The results revealed that there are appreciable 
differences in the statistics of temperature between SFT and CFT. 

The kinetic energy spectra follow the k~ 5 / 3 power law, and the Kolmogorov constant 
C K is 2.17 in SFT and 1.05 in CFT. The Helmholtz decomposition on the velocity showed 
that in SFT the compressive component of the kinetic kinetic energy spectrum in SFT 
obeys the k~ 5 / 3 power law, while in CFT it defers to the k~ 2 power law. The temperature 
has the Kolmogorov spectrum in SFT and the Burgers spectrum in CFT. The 2D contours 
of temperature and temperature dissipation in SFT display the sufficient mixing of the 
large-scale ramps and the small-scale cliffs. By contrast, the same 2D contours in CFT 
are dominated by the large-scale motions of rarefaction and compression. The major 
contribution to the power-law region of large negative dilatation are from preshocks and 
weak shocklets rather than strong shock waves. Our results showed that the power-law 
exponents are —2.5 for SFT and —3.5 for CFT. Using a theoretical model to handle 
the conditional averages of pressure and viscous, we obtained —2.56 for SFT and —3.43 
for CFT, which are close to the numerical values. The p.d.f. of temperature increment 
is concave and symmetric. When scale is comparable to the integral scale, in SFT it 
becomes Gaussian but in CFT it is still super-Gaussian. Moreover, unlike that in SFT, 
the rescaled p.d.f. in CFT collapses to the same distance, indicating that the scaling 
exponents for the statistical moments of temperature increment should saturate at high 
order numbers. 

We further studied the isentropic approximation in thermodynamic variables. Within 
a tolerance of 10%, in SFT the pointwise values of temperature and pressure is only 
1.6% failing to satisfy the isentropic relation. However, in CFT the deviation increases 
to 18.9%. For large negative dilatation, the ratio of T/(T) to (p/ (p)) 1 ^ 1 ^ 7 is near unity, 
meaning that the thermodynamic process is approximately isentropic in the compression 
region. Then, the accuracy of isentropic approximation was measured by introducing the 
quantity of p = p / (p) /(p/(p)) 7 . It showed that in SFT the p.d.f. of p is close to the 
Dirac delta function at M t = 0.1, and deviates from that as M t increases. At M t ss 0.6, 
p of CFT displays a broader distribution than that of SFT. Furthermore, it showed that 
the tails of the p.d.f. of p, Q(p), from SFT are well described by a theoretical model 
based on the Gaussian assumption of temperature distribution. 

The description in the statistical properties of temperature gradient revealed that the 
temperature tends to tangent to the vorticity in highly compressible turbulence. Com¬ 
pared to that in SFT, the temperature gradient in CFT is preferentially perpendicular 
to the solenoidal component and preferentially aligns with the compressive component 
of the anisotropic strain rate tensor. There are clear tendencies for the temperature gra¬ 
dient to align with the first eigenvector corresponding to the most negative eigenvalue, 
and to be perpendicular to the second and third eigenvectors. The conditional magni¬ 
tude average of temperature gradient has the teardrop shape, and is substantial in the 
region where the second and third invariants of the anisotropic velocity gradient tensor 
are respectively positive and negative. 

By employing a ” coarse-graining” approach, we studied the cascade of temperature 
in compressible turbulence. It was found that throughout scale ranges, the transport of 
temperature fluctuations is increased by the viscous dissipation at small scales and the 
pressure-dilatation at moderately large scales, however, is decreased by the SGS temper¬ 
ature flux, which preferentially occurs in the orientation where the temperature gradient 
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anti-align with the SGS velocity-temperature coupling. The appearance of plateau in the 
SGS temperature flux indicates the conservation of temperature cascade from large to 
small scales. The slope values for the cospectrunr of pressure-dilatation is —1.30 in SFT 
and —2.38 in CFT, indicating that the pressure-dilatation converges and be independent 
of scale at high enough wavenumbers. It provided the picture that the conversion between 
the kinetic and internal energy by the pressure-dilatation mainly occurs over moderately 
large scales of limited extent. The positive and negative components of pressure-dilatation 
are substantial at small scales. Once taking global averages, they basically cancel each 
other and make the outcome small. This does not contradict to the fact that the mo¬ 
tions of rarefaction and compression happen at all scales in compressible turbulence. The 
strongly positive skewness of the p.d.f. of pressure-dilatation from small scales implies 
that the conversion from kinetic to internal energy through compression is more intense 
than the inverse process through rarefaction. The 2D contours showed that in SFT the 
pressure-dilatation for large scales and the residual part from small scales distribute ran¬ 
domly, while in CFT they are concentrate in the rarefaction and compression regions. 
Finally, we observed that the fluctuation component of pressure-dilatation in SFT is well 
described by the Zernan model. 

The current investigation reveals a variety of unique statistical properties of tempera¬ 
ture in compressible turbulence, relative to the features of temperature in incompressible 
flow. These new findings can be largely understood through the effects of shock topology 
and the degree of compressibility. We limit our study to numerical simulation, a further 
work addressed the theoretical models for temperature in compressible turbulence will 
be performed in the near future. 
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